A Thermodynamic Model of Electric-Field-Induced Pattern 
Formation in Binary Dielectric Fluids 

M.D. Johnson,Q X. Duan, Brett Riley, Aniket Bhattacharya,[] and Weili Luo| 

^ ! Department of Physics, University of Central Florida, Orlando, FL 32816-2385 

> • (Dated: November 20, 2002) 

O ' 

Abstract 

£N| ' An electric-field-induced phase transition and pattern formation in a binary dielectric fluid layer 

are studied using a coarse-grained free energy functional. The electrostatic part of the free energy 
is a nonlinear functional of the dielectric function, which depends in turn on the local colloidal 
concentration. We determine the phase co-existence curve and find that beyond a critical electric 
field the system phase separates. Accompanying the phase separation are patterns similar to 
those observed in a spinodal decomposition of an ordinary binary fluid. The temporal evolution 

o. 

t O , of the phase separating patterns are discussed both analytically and numerically by integrating a 

Cahn-Hilliard type of equation. 
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I. INTRODUCTION 



Systems undergoing a phase separation from a disordered phase to a more ordered state 
often form interesting patterns as they select and approach the final state 0]. The process of 
phase ordering, or coarsening dynamics, can be classified into different universality classes 
depending on the physical dimension, the dimension of the order parameter, and conservation 
laws0, [|. Typically the associated scaling properties and dynamical universalities can be 
understood using a coarse-grained free energy functional. A time-dependent Ginzburg- 
Landau approach then allows study of the evolving patterns. 

Ferro- fluids H, which are colloidal suspensions of magnetic particles in solution, exhibit 
complex labyrinthine pattern formation when trapped between closely spaced glass plates 
and subjected to a transverse magnetic field||, |], |§. The complex patterns result from 
competition between the surface tension and long-range forces, and have much in common 
with patterns observed in amphiphilic systems and Langmuir-Blodget films. In this paper 
we report on a similar effect in colloidal suspensions in which phase separation and a con- 
comitant pattern formation are driven by an external electric field. Our results are based 
upon a macroscopic thermodynamic model of a dielectric binary liquid. This work was mo- 
tivated by recent experiments demonstrating such behavior 0, §j, and shares some features 
in common with the experiments. We reported some preliminary results in Ref. |9]. 

We consider colloids consisting of nanometer scale (typical diameter ~ 10 nm) dielectric 
particles suspended in a dielectric solution. In the present work we are interested in patterns 
on length scales much larger than the particle size and the typical inter-particle spacing. 
Consequently we use a coarse-grained description similar to the model B used to study 
phase separation in binary liquids^], [TOf] . For an order parameter we use the local particle 
concentration or volume fraction c(r) (0 < c < 1). We model the particle interaction using 
a local dielectric constant e(c(r)) which depends nonlinearly on the local concentration. 
The use of a phenomenological dielectric constant with a positive curvature is an important 
ingredient of the model. It turns out that beyond a critical value of the electric field the 
electrostatic energy dominates and drives phase separation and pattern formation. 

The approach developed in this paper is thermodynamic, based on a macroscopic free 
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energy. Analyzing equilibrium properties yields a phase diagram which exhibits instabilities 
that can lead to pattern formation. Some understanding of the dynamics of pattern for- 
mation then comes from investigating the corresponding time dependent Ginzburg-Landau 
equation. We begin by describing the model free energy and its equilibrium properties, and 
then turn to dynamics. 



II. MODEL FREE ENERGY 



Consider a binary dielectric liquid as described above placed between parallel capacitor 
plates separated by d. The liquid has capacitance 

C=-L [d\e(c(r)). (1) 



And 2 _ 

When the internal energy is dominated by the electrostatic contribution, the Helmholtz free 
energy is 

F = \QV - TS, (2) 

where V is the voltage between the plates, T the temperature, and S the entropy. The 
equilibrium state is obtained by minimizing F at constant Q, T. For our purposes it is more 
convenient to make a Legendre transformation to a different free energy: 

n = F -QV = -\QV -TS. (3) 

Equilibrium is found by minimizing Q at constant V, T (note that V is voltage and not 
volume). We will refer to Q as "the" free energy. 



For the dielectric constant e we use a phenomeno logical model due to LichtenekerP, ITT 



e(c) = £/e 7C , where 7 = hie p /ef. (4) 

This interpolates between the dielectric constant of pure solvent {ef) and that of pure particle 
(e p ) as the concentration varies from c = to c = 1. In this paper we use a physically 
reasonable value 7 = 30, ||. The specific choice of e(c) in Eq. (f|) is not essential. The 
results we report turn out to be largely independent of the form of e(c), as long as it has 
positive curvature (a point we discuss later). 
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We use an approximate entropy functional appropriate for interacting particles at low 
concentration |7|, |8|, [12|, [T3|j : 



S=~ J d 3 r [c(r) In c(r) + (1 - c(r)) ln(l - c(r))] . (5) 

This expression can be obtained by a simple counting argument. In the absence of interac- 
tions v is the particle volume and Eq. @ is the entropy of free particles |T3J]. Here however 
v is instead the so-called 'correlation volume,' which accounts phenomenologically for the 
particle interactions. 

In a typical experiment the plate separation d is so small that the concentration c(r) 
can be treated as effectively two-dimensional. Then we can let the coordinate r = (x, y) 
denote lateral position within the plane. In this limit the electric field E becomes uniform 
(E = V/d), and the free energy becomes 



Q = d J d r/(c(r)), where (6a) 
/( c ) = -^ £ (c) + —[c\nc+(l-c)\n(l-c)}. ( 6b ) 

07T V 

For e(c) given in Eq. (|3j), the temperature and electric field naturally scale to the dimen- 
sionless form 

f = kBT/v . (7) 

III. EQUILIBRIUM PHASE SEPARATION 

We are interested in potential phase separations in the system. These can occur if f(c) 
possesses inflection points|15[. Suppose that in equilibrium the system phase separates into 
volume fractions with concentrations C\^. If the average concentration is c, then these 
must satisfy 

2 2 

= i, ^2x^ = 0, (8) 

i=l i=l 

and so 

c 2 — c C — Ci 

xi = , x 2 = . (9) 

C2 — Ci C2 — Ci 
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To find the conditions under which the equilibrium state is phase-separated in this manner, 
we minimize the free energy £V Xif(ci) with respect to ci, c 2 , X\, x 2 , subject to the constraints 
Eq. (^|). Using Lagrange multipliers, this means minimizing 

2 

n = ^2x i \f{ Ci )-tic i -x]. (io) 

i=i 

Except at the endpoints (xi or q equal to or 1) this minimization yields 

f(ci) = /iCi + X (11a) 
f'( Ci ) = n. (lib) 

If ci and c 2 both satisfy Eq. ( |lla| ), then 

/(c 2 )-/(c 1 ) = /i(c 2 -c 1 ). (12) 

This result is summarized by a simple geometrical construction illustrated in Fig. |l|. If the 
straight line joining the points (ci,/(ci)) and (c 2 ,/(c 2 )) is tangent to /(c) at both points, 
then phase separation will occur for c\ < c < c 2 . This can occur only if f(c) has two or 
more inflection points (see Fig. |l|). For our model, the endpoints c = 0, 1 never correspond 
to a minimum free-energy state. 

Consider now the free energy density f(c) for our dielectric fluid, Eq. (|6l5|). The entropic 
term —TS has positive curvature. At high temperatures this term dominates, f(c) has 
positive curvature everywhere, and the equilibrium state is homogeneous. This is illustrated 
by the highest-temperature trace (T = 20) in Fig. |2|, which uses the model dielectric function 
Eq. (|4]). As the temperature is lowered the electrostatic energy term —QV/2 = —E 2 e(c)/8ir 
plays an increasingly important role. When e(c) has positive curvature, the electrostatic 
term has negative curvature, and at low temperatures this term leads to inflection points 
and phase separation. This is illustrated by the two lower-temperature traces in Fig. |2|. 

The phase diagram for the free energy Eq. (|6bD , developed using the above procedure, 
is shown in Fig. |^ for 7 = 3. Below the coexistence curve (the solid line) the homogeneous 
phase becomes metastable and in equilibrium the system phase-separates. This curve is the 
locus of points obtained from solving Eq. (|TT| ) at all temperatures. Although there is no 
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exact expression for the coexistence temperature, for the physically important case of low 
concentrations we find approximately 

~ _ (l + 7- 7 c)e^ c -e^ 
Inc 

For 7 = 3 this is accurate for concentrations c < 0.3. 

Lying below the dashed phase boundary in Fig. |3] is the classical spinodal region, where 
the homogeneous state becomes unstable [/"(c) < 0]. This phase boundary is given by the 
condition /"(c) = 0, which yields 

T s = 7 V c c(l-c). (14) 

The coexistence and spinodal curves touch at a critical point (marked by a diamond in 
Fig. |3|). This occurs at the critical concentration c c and critical temperature T c where the 
two inflection points present at lower temperatures coincide [and so f"(c c ) = = f"'(c c )]. 
For the model dielectric function Eq. (|j) the critical values are: 

c c = (7 - 2 + V7 2 + 4) (15a) 
T c = 7 V Cc c c (l-c c ). (15b) 

The instabilities give the possibility of pattern formation, which we now discuss. 



IV. PATTERN FORMATION: DYNAMICS 



Pattern formation after quenches into the metastable or spinodal regions provides impor- 
tant insight into the underlying mechanisms for phase transformation in the system. For a 
system like the one modeled here it is experimentally easier to increase the voltage suddenly 
while holding the temperature fixed; by Eq. (0) this is equivalent to a quench. Here we 
investigate pattern formation under quenches using a phenomenological diffusion equation 
developed from the free energy. 

Notice first that when inhomogeneities arise the free energy needs an additional |Vc| 2 
term. This can be viewed as the surface energy between domains of different concentration; 
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such terms arise generally in any interacting system. Consequently in this section we write 
the free energy 



d 3 r 



/(c(r)) + i« | Vc| 



(16) 



where k is a phenomenological parameter. 

Fick's law relates the concentration current density to the local chemical potential fi: 



J r 



(17) 



where M is the mobility (taken to be constant) and vq is the volume of the colloidal particle. 
Combining this with the continuity equation c + V • j = yields a diffusion equation 



dc 

— = Mv V 2 fi. 

The local chemical potential is obtained from the free energy: 

5Q 



(18) 



M r ) = 



5c(t) 



v [/'(c(r)) - ^V 2 c(r)] . 



(19) 



TE 



Inserting this into Eq. (|TJ) gives 



dc 



- = M^[V 2 f(c(r))- K V 4 c(r)]. 



(20) 



The free energy density f(c) is in Eqs. (|6b"|,|4]). A convenient dimensionless form comes from 
writing time in units of r and length in units of A, where 



A 



Then the diffusion equation becomes 



K ( 8tt Y 
8ttk 



(21a) 
(21b) 



dc 
Of 



^ = V 2 



/'(c) - V 2 c 



where 



/(c) = -e 7C + f [c In c + (1 - c) ln(l - c)] . 



(22a) 



(22b) 



7 



Eq. ( |22a| ) is nothing but the time- dependent Ginzburg-Landau equation for a scalar con- 
served order parameter; however, the free energy in this particular case is constructed with 
the parameters directly relevant to experiments on dielectric colloidal fluids [0, §]. We study 
pattern formation after quenches by solving this nonlinear diffusion equation, starting from 
a uniform distribution (plus a small random component). 

One can get some intuition about dynamics by linearizing the diffusion equation^]. The 
result is accurate for early times, before the nonlinear terms become significant. Write 
c(r) = c + Sc(r), where c is the average concentration. For small 6c expand 

/(c + 6c) = /(c) + faSc - l -b (5c) 2 + ..., (23) 

where 

fa = f(c) = - 7 e 7C + f [lnc - ln(l - c)] (24a) 

b = -f"(c) = (24b) 
c(l - c) 

T s is given in Eq. ([141) . The linear and constant terms in Eq. (|23"D play no role in the 
dynamics. Keeping only the second order term from Eq. (|23|), the free energy Eq. (|T6"D 
becomes 

n^^Jd 3 r [k\\75c\ 2 - b(5c) 2 ] . (25) 

This expansion connects our model to the coarsening dynamics of a binary alloy described 
by a conserved scalar order parameter 0(r). The free energy for such a system is generally 
assumed to have the Ginzburg-Landau form 

F= 1 -Jd 3 r [Km 2 - b<J> 2 + U -<f\ (26) 

where k, and u are assumed to be positive. Below the critical temperature b becomes 
positive, yielding a broken symmetry. Higher-order terms beyond 4 are irrelevant for the 
dynamic universality class. In fact, the early-time dynamics are dominated by <fi 2 terms. 
From the correspondence between Eqs. (p5[), (f26|) we expect our model to have early-time 
dynamics very similar to those of model B. 



Using the expansion Eq. fl25|) linearizes the diffusion equation Eq. (|22|) 

d{5c) 



-V 2 [b5c + VHc]. (27) 
This is solved by expanding the fluctuation 5c in Fourier components: 

Sc{*,t) = J (0 elk Fc (k,t). (28) 

We find 

c(k, t) = c(k, 0)e akt , where a fe = k 2 (b- k 2 ). (29) 

Outside of the spinodal region T > T s and so by Eq. ( |24b| ) 6 < 0. Consequently all 
components c(k, t) decay. This describes the homogeneous equilibrium state: fluctuations 
away from a uniform concentration decay. 

However, within the classical spinodal region b > 0. Thus long wavelength modes, those 
with wave vector k < Vo, grow exponentially. The most rapid growth occurs for k m = y/b/2, 
where reaches its maximum value b 2 /4. Thus in the classical spinodal regime, where the 
homogeneous state is unstable, early-time exponential growth leads to structures a typical 
(dimensionless) length scale A;" 1 ~ (T s — T) -1 / 2 . These structures grow exponentially in 
time with a (dimensionless) time scale l/«fc m = 4/6 2 ~ (T s — T)~ 2 . 

For example, consider the structure factor 

S{k,t) = (c{k,t)c{-k,t)). (30) 

Here angle brackets represent a statistical average. Suppose that the initial fluctuations 
c(k, 0) are small and uncorrelated: 

(c(k,0)c(k',0)) = A 2 ( 5 2 (q-q'). (31) 

Then for early times 

S(k,t) = A 2 e 2akt . (32) 

This grows ever more peaked at k m as time increases. Fourier transforming gives the real- 
space correlation function 

(Sc(v,t)5c(r',t)) ~ -^e b2 '/ 4 J (Ur-r'|) (33) 



which exhibits the length and time scales described above. 

Including dimensions, the structures developing in early times have a length scale 



km 




T s -T 



and grow with a time scale 



4r 



K ( 8tt V 
Mv$ \e f E 2 ) 



2c(l 



Mv 2 



2c(l - c)k B /v 
T-T s 



(34) 



(35) 



T-T s 

Deep quenches into the classical spinodal region should lead to dynamically-developing 
structures on a length scale L given by Eq. fl34l) . These structures result from a long- 
wavelength instability. They represent not the equilibrium state, but the early-time evolution 
from the highly non-equilibrium homogeneous state toward the (eventual) phase-separated 
equilibrium. The linearized description is accurate for very small initial times. Soon the fluc- 
tuations grow large enough for the nonlinear terms to play a role, and we turn to numerical 
means for the solution of Eq. (|22|) . 



V. NUMERICAL INTEGRATION OF THE DIFFUSION EQUATION 



We numerically integrate the lattice discretized version of Eqs. (|22|) using a first-order 
Euler scheme, choosing step sizes to avoid unphysical instabilities^]. This algorithm is 



adequate for times during which the system begins to organize into patterns. To carry 
the simulation all the way to phase separation requires a more sophisticated approach ||i~7 



However the model we use also omits hydrodynamic effects which may be important for the 
late time coarsening process. Consequently here we present results during the era of pattern 
formation only. 

The initial configuration consists of a uniform concentration c plus a small random com- 
ponent. Starting from this initial state, Eq. (^) is integrated for various values of 7 and T. 
In agreement with the results of the previous sections we find that phase separation occurs 
below the coexistence temperature T coex (or, equivalently, above a threshold electric field). 
In the phase-separation region the concentration develops labyrinthine patterns consisting 
of stripes of higher and lower concentration. 
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Typical results for a quench deep into the phase separation region are shown in Figs. |] 
and [5[ The latter clearly shows that a labyrinthine pattern has developed. 

We monitor two quantities during the simulation. At each time step we search for the 
maximum and minimum concentrations, and define a contrast parameter (c max — c m i n )/2c. 
This parameter indicates the visibility of the developing pattern. Another quantity of inter- 
est is the pair correlation function 

G(r,t) = (hj W,t) - c)(c(f+ f,t) - c)d 2 ?\. (36) 

The average domain width R g can be calculated from the first zero of the normalized pair 
correlation function G(f, t) = G(f,t)/G(0,t). The upper panel of Fig. |4] shows the early time 
variation of the contrast parameter and the domain width R g as a function of the scaled 
time. The average domain width after a rapid initial increase saturates at this temperature 
to a value ~ 4.5A. 



VI. SUMMARY 



We have studied the electric-field-induced phase separation of a dielectric binary fluid 
using a coarse-grained free energy functional. The free energy investigated is guided by 
experiments on field-driven instabilities in ferro-colloids, and shares features of a Ginzburg- 
Landau free energy that describes the phase separation of a binary mixture. Certain other 
features of ferro-particles which would require a more complicated order parameter have 
been omitted, e.g. the orientational ordering arising from the magnetic dipolar interaction. 

The first order co-existence curve of the model clearly shows how a spinodal decomposition 
can be driven by an electric field in colloidal systems. Analytical estimates of the critical 
field are consistent with the numerical results. Deep quenches show a typical labyrinthine 
pattern formation with characteristic length and time scales set by the free energy. 

The phenomenology used here is most accurate at low particle concentrations. In par- 
ticular, the entropy is essentially that of free particles, modified to take interactions into 
account via the phenomeno logical correlation volume v. The resulting entropy is at least 
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qualitatively correct at low concentrations but not high. Experiments with low overall con- 
centrations should be described with reasonably accuracy by this model. 
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FIG. 1: Two points on a curve can have a common tangent only if the curve has at least two 
inflection points. The inflection points at ca,cb are indicated by filled circles and the phase 
separation points at c\,C2 by open circles. 
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FIG. 2: (a) f(c) vs. c for the three scaled temperatures T = 1, 10, 20, using Eqs. (6b, |) with 7 = 3. 
At the two lower temperatures two inflection points are present. These are more easily seen in (b) 
f"(c) vs. c. The inflection points are indicated by filled circles. 
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FIG. 3: Phase diagram for the free energy given in Eqs. (|6b|j4D for 7 = 3. Temperature is T. Shown 
are the coexistence curve icoex (solid) and the spinodal T s (dashed). The critical point Eq. (|15| ) is 
marked by a diamond. 
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FIG. 4: Early stage increase in contrast (c max — c m i n )/2c and domain size (inset) as a function of 
time. Here time and space are in the units of Eq. (|2l]). The simulation was done on a 128x128 
lattice with an initial random configuration ±0.01 around c = 0.1 for 7 = 1.0 and T = 0.045. 
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FIG. 5: A typical snapshot at time 20. Concentrations greater (less) than the average c are dark 
(light). 
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